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ABSTRACT 

This paper presents recent time-dependent measurements of neutral B meson 
oscillations. Similar to the K'^-K^ system, there are two such systems involving 
the b quark: B^~B^ and B^-B^. Thus the physical states are respectively 
Ks and Kl, {Bd)s and {Bd)L, and {Bs)s and {Bs)l- The oscillation between 
each pair of states can be used to determine their mass difference. The present 
world average for the {Bd)s~{Bd)L mass difference is Anid = 0.457±0.019 ps""'^ 
(or (3.01±0.13)xl0-'' el/). Using fss = 12% (the fraction of B^ produced in 
b events), the current lower limit on the corresponding Amg is 6.1 ps~^ (or 
4.0x10-3 ey). 

1. Introduction 

Since there are three known quarks of charge — ^, namely d, s and b, there are 
three similar neutral particle-antiparticle systems: 

K^{sd)-K%sd), Bd{bd)-Bd{bd), and B,{bs)-B,{bs) . 

Of these three, the K^-K^ system is best measured experiment ally0i and understood 
theoretically!. However, the theoretical analysis applies equally well to all three cases. 

As the topic here is the BB mixing, let B and B denote the flavor state in all 
three cases, i.e., 

B = Bd or 5s 

and B = Bd or B,, 

while the corresponding weak eigenstates are 

^5 = Ks, iBd)s or {B,)s 
and Bl = Kl, {Bd)L or {B^)l. 

Because of CP non-conservation, Bg and B^, which are not orthogonal, can most 
generally be related to B and B by 

Bs = {\v? + \q?y~' {vB + qB) 
and Bl = {\pf + Igf)'^' {pB - qB) . 

Let r = (Ts + Tl) /2 and m = {ms + trl) /2 be the average width and mass of 
Bs and Bl, while AF and Am are the differences 

Ar = - > 
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and 

Am = \ms — 171^1 . 

Let VB,uif) and V^^ui't) be the probability distributions for a meson which is created 
as B (or B) to decay as a 5 (or B) after a proper time t, and VB,m{t) and VB,m{^) be 
those for a meson created as B (or 5) to decay as a S (or 5), where the subscripts u 
and m stand for unmixed and mixed respectively. These four quantities are given by 

|2 



Vb,u (t) 



VB,m{t) 



VB^uit) 



VB,m{t) 



+ 



r2-(Ar/2)'^ ' r2+(Am) 

m 



-rt 



+ 



ipi -kr 



r2-(Ar/2)'' ' r2+{Am) 







2 


r 


r ipp+kP 




|2 1 |2 n 

\p\ -kl 




Lr2-(Ar/2)2 


r2+{Am)2j 



-r* 



r2-(Ar/2)2 



Ar 

cosh 1 + cos Am t 



cosh 1 — cos Am t 

2 



cosh t + cos Am t 



cosh t — cos Am t 

2 



r2+{Am)2_ 

Although these three systems, K^-K^, B^-Ba, and Bg-Bg can be described by the 
same set of equations, the different decay modes lead to significant differences in the 
behavior of these systems. For K^-K^, the vrvr mode dominates with the immediate 
consequence that Ts ^ Tl, in fact. 
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In contrast, both Ba and Bs have many important decay modes. Indeed, for both 
cases the width difference AF comes from decay modes that are available to both B 
and B. Using the many measured branching ratios for B^, a generous estimate gives 
(AF/F)^ < 5%, perhaps much less. Since there is very little experimental information 
about the decay modes of Bs, no such firm statement can be made about the ratio 
(AF/F)g, but it is also believed to be small. In this talk, the width difference AF will 
be neglected for both the Bd-Bd and the B^-Bs systems. 

From Eq. (|1]), if the effect of CP non- conservation is neglected, then the formulas 
for the four probabilities simplify to 
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If, furthermore, AF is neglected, the preceeding equations can be written as 

F 
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e~^* [1 + cos Ami] 
e~^^ [1 — cos Ami] 
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For the K^-K^ system, the mass difference Am was measured a long time ago! to 
be 3.51 X IQ^^eV . The purpose of this talk is to present and discuss recent experimen- 
tal results on Am^ and Amg for the Bd-Bd and Bg-Bg systems. The most sensitive 
measurements of Am^ and Amg are obtained through time- dependent measurements, 
which investigate Vu and Vm directly, and this talk will consider only time-dependent 
measurements of these quantities. Since neither CP non-conservation nor the width 
differences AF have been observed in these systems, the analyses have been carried 
out on the basis of Eq. @ rather than the more accurate Eqs. (|l|) or (|^). 



2. Experimental Overview 

In order to perform a time- dependent measurement of Am^ or Amg, one must 
measure the proper time of the decay of the B meson, and determine its production 
flavor (i.e., _B or 5 at production) and decay flavor (i.e., _B or _B at decay) in order 
to ascertain whether the B meson is mixed or unmixed. The value of Am is then 
found from the fraction of events identified (or tagged) as mixed or unmixed as a 
function of the measured proper time. Fig. |l] addresses the experimental sensitivity 
of such a measurement. Fig. [l|(a), (b) and (c) each show the decay probabilities 
Vu(t) and Vm(t), discussed above. Each figure assumes a B lifetime of 1.5 ps, and 
shows the effect of different values of Am on Vu(t) and Vm(t). Fig. |l|(a) shows these 
probabilities for Am^^ = 0.5 ps~^. It demonstrates that, since the typical experimental 
resolution for the LEP experiments is about 0.25 ps in B meson proper time, it is 
a relatively easy job to measure Am^ due to the large oscillation period. Fig. |I|(b) 
illustrates that if Arrig = 5 ps~^, the oscillation period is still within a comfortable 
reach of the experimental sensitivity. For Amg = 15 ps~^. Fig. |l](c) shows that 
the experimental sensitivity for LEP experiments makes it difficult to extract the 
frequency of oscillation. 

The proper time t of the B decay is obtained through 

where / is the decay flight distance between the B production point and decay point, 
and rriB and pb are respectively the mass and momentum of the B meson. The 
flight distance is measured with the aid of silicon microvertex detectors which allow 
the production and decay vertices to be reconstructed precisely. The B meson decay 
length is determined by reconstructing a decay vertex formed from a lepton with high 
transverse momentum (or pt) and a "charm track". The "charm track" is formed by 
combining information from several tracks which are not consistent with coming from 
the production point of the B and form a secondary "charm vertex". (In the case 
where a D* from B decay is fully reconstructed, a variation of this method is used; see 
section 3.3). Because of the presence of tails in the flight distance resolution, it must 
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Fig. 1. Vu and Vm for (a) Am = 0.5 ps"^. (b) Am — 5 ps~^. (c) Am ~ 15 ps~^. The solid line 
shows the exponential decay of the B meson with a lifetime of 1.5 ps. The dashed line shows the 
Vu distribution, and the dotted line shows the Vm distribution. 



be parametrized with several Gaussians. Typically half of the measurements fall in 
the "core", where the error is smallest. This core resolution is 260 iim for alephS, 
340 Iim for DELPni, and 400 /im for OPAll. 

The B momentum is obtained by reconstructing the momenta of the charged 
and neutral decay products of the B. The charged momentum can be reconstructed 
by simply summing the momenta of charged tracks consistent with coming from the 
decay of the B. These usually include a lepton and other charged tracks from a charm 
meson decay vertex. The neutral energy reconstruction is generally more complicated, 
involving information from the whole event, the beam energy, and energy- momentum 
conservation. The B momentum core resolution in ALEPH0, DELPHli and OPAlI 
is about 8-10%. 

The charge of the h quark when it is created (the production flavor) is typically 
found in one of two ways. Some analyses require a lepton in the hemisphere oppo- 
site to the lepton used to determine the decay flavor, and use its sign to determine 
the production flavor. Others use Jet charge techniques, which weight momentum 
information from charged tracks in the event to determine the production flavor. 
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The charge of the b quark when it decays (the decay flavor) can be measured in 
a variety of ways. In some analyses, the sign of a high pt lepton is used to identify 
the decay flavor. In other analyses, a D*^ is reconstructed, and the sign of the D*^ 
is used to determine the decay flavor. 

Details of the speciflc methods used are discussed in the next section. 



3. Measurement of Am^. 

Bd oscillation was flrst observed by alephI at LEP in 1993. Since then, a vast 
number of measurements have emerged using a variety of methods. The popular 
methods are described here, namely the "Lepton- Jet charge" method, the "Lepton- 
Lepton" method and methods using a D*. The names of these methods are chosen 
such that the word before the hyphen refers to the way the decay flavor is determined, 
while that after the hyphen is for the corresponding production flavor. 



3.1. Lepton- Jet charge Method. 

In the Lepton- Jet charge method, events with semileptonic decay b —>■ Xi~u or 
b — > Xi'^u, {i"^ = e^or^M^) are selected with a high pt lepton; the charge of the lepton 
from these decays identifles the decay flavor of the b quark. Leptons from other 
sources, particularly cascade decays 6 — >• c — dilute the sample of b ^ £~ , but the 
high Pt lepton is nevertheless a good measure of the b decay flavor. The production 
flavor is then tagged by a Jet charge technique, discussed below. The B meson decay 
vertex is determined by the intersection of the high pt lepton and the "charm track" 
on one side of the event, as described in section 2. The schematic of this method is 
shown in Fig. |](a). 

There are several different algorithms used to determine the b production flavor. 
The choice of charged tracks used in these analyses varies. ALEPII0 and DELPHIH use 
only the tracks in the hemisphere opposite to the lepton, while opal! uses charged 
tracks from both the lepton jet, and the highest energy jet which does not contain the 
lepton, also called the opposite side jet, in calculating its jet charge. For convenience, 
this paper will refer to both jet charge and hemisphere charge measurements as jet 
charge measurements. 

The weighting scheme for the ALEPH and DELPHI results takes the form 



i=l WiQi , , 



where uh is the number of tracks in the opposite hemisphere, is the charge of the 
track, and Wi is the weight used, taking the form \pi ■ e|'^, where e is the direction 
of the thrust axis for ALEPH, and the direction of the sphericity axis for DELPHI. 
ALEPH calculates this weight using k = 0.5, while DELPHI uses k = 0.6. OPAL 
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uses a different jet charge, defined by 



g2j = E^«-ioE?n^) (6) 



J J -^beam 

where the first sum is over the tracks in the jet containing the high pt lepton, and the 
second sum is over the opposite side jet. In Eq. (|^), -Ebeam is the beam energy, and 
Pi\\ is the charged track's momentum parallel to its jet axis. 

Because only one high pt lepton is required, the "Lepton- Jet charge" method 
retains a relatively large sample of events and hence possesses a strong statistical 
power. The tag rate, or the fraction of events correctly identified as mixed or unmixed, 
for jet charge analyses is approximately 70% for both mixed and unmixed events. The 
ALEPH result studies the time dependence of the lepton-signed jet charge (jet charge 
multiplied by sign of lepton) distribution, without explicitly identifying events as 
mixed or unmixed. The results of Am^ from alephI, Delphi! and opal! usmg 
this method are shown in Figs. |^(b), |^(c) and H(d). Throughout this report, where 
errors are given, the first is statistical and the second is systematic. 



3.2. Lepton-Lepton Method 

Dilepton measurements are in many ways similar to jet charge measurements. In 
the Lepton-Lepton method, events with semileptonic decay b Xi~u or 6 — > Xi'^u, 
{i^ = or /i^), on both sides of the event are selected. The B meson decay vertex 
is determined by the intersection of the high pt lepton and the "charm track" on one 
side of the event, as described in section 2. The B meson decay flavor is tagged by 
the sign of the high pt lepton on the flight distance side of the event, just as in the 
Lepton- Jet charge method. The flavor at production time is tagged by a lepton in 
the opposite hemisphere. Contributions from mixing of the opposite side B hadron 
are independent of the proper time in the flight distance hemisphere, and their effect 
is factored into the tag rate calculation. The entire process can then be repeated 
with the roles of the leptons reversed, giving up to two measurements per event. The 
schematic of this method is shown in Fig. ^(a). 

Because of the requirement that a lepton be found in each hemisphere, Lepton- 
Lepton measurements have a smaller event sample than Lepton- Jet charge analyses. 
Compensating for their smaller event sample, Lepton-Lepton analyses have superior 
tag rates, correctly identifying events as mixed or unmixed 80% of the time. This 
gives them sensitivity comparable to the Lepton- Jet charge method. The results from 
ALEPhIzI, OPAlI, and CDF0 using this method are presented in Fig. Kb), ||(c), and 
1(d). 

Using the Lepton-Lepton method, the DELPHI experiment^ finds Am^ = 0.42 ± 
0.08]'^o;o7 ps~^. DELPHI extends this method by including the use of a charged kaon 
to identify the decay flavor of the B meson, making use of DELPHI'S unique feature. 
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Fig. 2. Measurement of Am^ with the Lepton-Jet charge method, (a) Schematic, (b) ALEPH 
measurement, (c) DELPHI measurement, (d) OPAL measurement. 
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Fig. 3. Measurement of Am^ with the Lepton-Lepton method, (a) Schematic, (b) ALEPH mea- 
surement, (c) OPAL measurement, (d) CDF measurement. 
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the RICH counters. Thus, in the flight distance hemisphere, the measurement uses a 
lepton or a charged kaon coming from the secondary vertex, relying on the dominant 
decay chain b ^ c ^ s to identify the B flavor. Such a kaon can be used in either the 
flight distance hemisphere to determine the decay flavor, or the opposite hemisphere, 
to determine the production flavor. This analysis also incorporates the jet charge or 
lepton in the opposite hemisphere to determine the production flavor. The DELPHI 
measurement using this Lepton-Kaon-Jet charge methodi is Am^ = (0.563to;o46 ± 
0.058) ps~^. Because of the strong statistical correlation expected between this result 
and the other inclusive DELPHI results, it has not been included in the final Am^ 
average given in Section 3.4. 



3.3. Methods using a D*. 

It is also possible to carry out time- dependent measurements of Am^ by recon- 
structing a D*^ from the decay of a -B meson. By selecting a charged D* sample, 
it is possible to obtain a very pure sample of Bd mesons. Though some B^ mesons 
contribute to the D*~ sample, the contamination is small, and the effect of these 
B~^ decays can be included in the fit for Am^. Because a D* candidate must be 
reconstructed, these analyses typically have much smaller event samples than either 
the Lepton- Jet charge or the Lepton-Lepton measurements. Two methods are de- 
scribed here, the "D*-Lepton or Jet charge" method, and the "D* Lepton- Jet charge" 
method. The schematic for these methods is shown in Fig. |^(a). 

In the D*-Lepton or Jet charge analyses, a Z)" sample is reconstructed using the 
decays to Ktc, Kinr^, and Kirmi. and then the D'^ is combined with a pion to produce 
a charged D*. The decay flavor of the B^ is identified by the sign of this D*. The 
production flavor can be identified with either a jet charge technique or with a lepton 
in the hemisphere opposite the D*, as discussed in the previous sections. As the 
pion from the D* decay has low momentum, and travels along the flight direction, 
it cannot be used to determine the B decay point. The apparent D° decay vertex is 
used to infer the B flight distance. The result from the ALEPH experiment^ using the 
D*-Lepton or Jet charge method is shown in Fig. ^(b). The DELPHI experiment^ 
obtains Am^ = 0.470 ± 0.086 ± 0.061 ps"^ with this method, while OPAL0 finds 
Arrid = 0.57 ± 0.11±0.02 ps~^ 

In the D* Lepton- Jet charge analyses, a D* sample is produced as described above, 
and a lepton in the same hemisphere is used to form a Bd decay vertex, for measuring 
the flight distance. The decay flavor is determined by the sign of the D*, and the 
production flavor is determined by a jet charge technique, as discussed in previous 
sections. The result from the DELPHI D* Lepton- Jet charge analysiJl is shown in 
Fig. 1(c), while the OPAL result§ is shown in Fig. §(d). DELPHI has performed 
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Fig. 4. Measurement of Arud with the D* methods, (a) Schematic, (b) ALEPH L)*-Lepton or Jet 
charge measurement, (c) DELPHI D* Lepton-Jet charge measurement, (d) OPAL D* Lepton- Jet 
charge measurement. 

The combined D* based resuh from DELPHI is Ama = 0.421 ± 0.064 ± 0.042 ps'^ 
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an average of their D* based analyses!, giving Am^ = (0.421 ± 0.064 ± 0.042) ps ^, 
which has been included in the LEP and world averages computed in Section 3.4. 

3.4- Average of Am a Results. 

The process of performing an average on the measurements of Am^ is complicated 
by the presence of correlations between the various measurements. These correlations 
can be statistical, coming from overlapping data samples, or systematic, coming from 
common assumptions used in making the measurement. 

Where measurements are statistically correlated, the degree of correlation is gen- 
erally unknown. To minimize correlations, when there have been multiple measure- 
ments of similar types performed on the same data sample, the average includes only 
those results with the smallest errors. Where measurements are systematically cor- 
related, it is possible to judge the degree of correlation between different results by 
looking at common correlated systematic errors. 

The main correlated systematic errors come from the lifetimes and fractions of 
individual B hadron species. Thus, the averaging process considers correlations re- 
lated to the lifetimes and fractions only. There are other errors which are correlated, 
in theory, such as the parametrization of the b fragmentation function, and the de- 
cay length resolution in individual LEP experiments, but these errors are generally 
smaller, and their correlations can safely be neglected. 

Taking these correlations into account correctly is complicated by the fact that 
each measurement parametrizes these systematic effects in a different way, and use 
different central values and errors on their parameters. The average considers the 
individual parametrizations of these different experiments, and performs a constrained 
for Arrid, the lifetimes and the B hadron fractions. 

Using this averaging technique, and excluding the ALEPH lepton-lepton mea- 
surement and the DELPHI Lepton-Kaon-Jet charge measurement due to statistical 
overlap, the LEP average is found to be Arrid = 0.458 ± 0.020 ps~^. Including results 
from CDF and time integrated measurements from T(4s) in the average yields Arria 
= 0.457 ± 0.019 ps~^. A summary of the results is presented in Fig. ^. 

4. Measurement of Arris. 

Both Am^ and Arris are mass differences between particles and hence they are of 
direct physical importance. A further motivation for measuring these quantities comes 
from the fact that these mass differences are due to high-order weak interactions. 
The important diagrams for these interactions are shown in Fig. ^ together with 
similar diagrams with the W and top quark lines exchanged. Computation from 
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Fig. 5. Summary of results for Am^. 
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calculation of Am^ and Amg. 
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where fisg and fjE^ are the QCD correction factors for the Bg and Bd, expected to 
be similar, and is the ratio of hadronic matrix elements for the Bg and B^. Es- 
timates bom lattice QCD0 a„d QCD sum rulefl are eonsistent with a valu-S of 
= 1-16 ± 0.10. Measurements of Arrid and Amg can therefore be used to determine 
|Vts/Vtd|. This ratio of the CKM matrix elements is of special interest because it 
appears in one of the most useful unitarity triangles given by 

^ + C + ^ = 0. (8) 

For a number of reasons, the measurement of Arris is much more difficult than that 
for Arud- The theoretical expectation that Amg is significantly larger than Am^ has 
been confirmed by experimentsi'0. As illustrated in Fig. |I], such larger values of Arris 
lead to rapid oscillation, which complicates the measurement. A second difficulty 
comes from the fact that /bs, the fraction of Bg produced in b decays, is substantially 
smaller than fs^ and is not well measured. 

Three methods for determination of Arris will be described here. They are the 
"Lepton-Jet charge" method, the "Lepton-Lepton" method and the "Lepton-Kaon 
Correlation" method. 



4-1- Lepton-Jet charge Method 

This method has been used by the ALEPhI, DELPHli, and OPAlI collaborations 
at LEP, and its schematic is that of Fig. ^(a) with Bd and Bd replaced by Bs and 
Bg. Again, the B decay vertex is formed by the secondary vertex including a high pt 
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Fig. 7. Lepton signed same side minus opposite side jet charge distributions for a) mixed B'^ mesons, 
and b) unmixed b hadrons in the full Monte Carlo simulation. Plotted in c) is the sum of all Monte 
Carlo events, normalized to the number of events in the data. The data points are superimposed 
with error bars. The vertical lines indicate the selection requirement \Qs — Qo\ > 0.2. 



lepton. The decay flavor is tagged by the high pt lepton while the production flavor 
is tagged by the jet charge technique. 

Since the lower bound from alephB using this method is the best one for Arris, it 
will be described in detail here. The main difference between the ALEPH Lepton- Jet 
charge method for ArUg and that for Am^ described in Sec. 3.1 is that a different jet 
charge algorithm is used. Instead of using the Qj^ in Eq. (^) where only the charged 
tracks in the opposite hemisphere are used, the new jet charge algorithm makes use 
of information from both the opposite jet and the flight distance jet. The weight 
applied to each track in computing the jet charge value for the event is the track's 
rapidity with respect to the jet axis. More precisely, define 

where S and O indicate the sum is over tracks in the same side jet (the jet with the 
high Pt lepton) and opposite side jet, respectively. The rapidity, Ui, is given by 

1 , Ei + Pill , , 
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Fig. 8. The tagged mixed fraction of events as a function of measured proper time, for a) Monte 
Carlo with Arris = 1-6 ps^^, and b) data events. The superimposed curve for the Monte Carlo in 
a) is the expected distribution for Arris — 1-6 ps~^. The solid curve for the data in b) assumes 
Amg = 30 ps~^, while the dashed curve in the insert is the expected distribution for Amg = 6 ps~^. 
The small proper time region of the plot is expanded to emphasize the part most sensitive to Bg 
oscillations. 



The jet charge variable used to identify the production flavor is then 

Q = Qs~Qo. (11) 

For the purpose of tagging mixed or unmixed events, the charge qe of the high pt 
lepton and the sign of the above Q are used. Events are identified as follows: 

H — , — h < — > mixed events 
+ + , < — > unmixed events 

As shown in Fig. a cut requiring |g£ ■ Ql > 0.2 is imposed. With this cut, the tag 
rate for unmixed events, Au, is about 80% and for mixed events. Am, is about 60%. 
Since there are eight times more unmixed events than mixed events, it is essential to 
have a high tag rate for unmixed events. This is a great advantage of using the Q 
defined in Eq. ([TT|). To make optimal use of the experimental data, the method of 
maximum likelihood is used to extract the value of Arris. 

Fig. ^ shows the tagged mixed fraction for Monte Carlo and data. The value 
of Arrid is determined as a check of the tag rates used in this analysis, and of the 



15 



12 

10 



M 8 
O 



o 



ALEPH (91-94) 



Phys.lett. B356 (1995) 409 
CERN PPE/95-84 



r a) \V 


fBs=12% 




— Data 

MC Average (Am=30ps"') 




• 95% CL points 


r • 









Oh 



6 



<1 

u 

in 

ON 



2 4 6 8 10 12 14 

Am (ps"^) 



b) 




6.1 


6.3 




6.5 ps'^ 




5.6 ,,- 








5.2 




















" , , , 1 




1 







10 12 14 



16 



18 



Fig. 9. Superimposed on the data likelihood curve are the average of the fast Monte Carlo AlogX 
values, and the 95% confidence limit points. AlogL is defined as the (— logL) value at the Amg 
value, minus the minimum (— log-L) value. A limit curve is drawn through the 95% confidence limit 
points. The data curve crosses the limit curve at Amg = 6.1 ps~^ for fraction /sg = 12%. Also 
shown is the average of the A log L curves for 200 fast Monte Carlo samples in which mixing is 
near-maximal (Amg = 30 ps~^). b) the results of 95% confidence level limit in Amg as a function 
of the fraction /sg- The limits for Amg are 5.2, 5.6, 6.1, 6.3, and 6.5 for = 8%, 10%, 12%, 
14% and 16%, respectively. 
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peformance of the likelihood fit. Assuming a B meson lifetime, r^, of 1.5 ps, 
fraction, /b^, of 0.4, B^ fraction, /^s, of 0.12, Amg of 30 ps^^, and Am of 0.6, a 
two-dimensional fit with the ALEPH data is performed for Am^^ and A^. This gives 
Au = 0.792 ± 0.003 and Am^ = 0.47 ± 0.04 ps~^, the latter value being in agreement 
with the world average of Fig. Similar fits with Monte Carlo give Au = 0.792±0.003 
and Anid = 0.48 ±0.05 ps~^ to be compared with the input values of 0.790 and 0.467 
ps^^ respectively. 

The ALEPH result is shown in Fig. ^. This figure shows the AlogL curve for 
the data as a function of Ams, where A logL is defined as the negative log likelihood 
value (— logL) at a given Am-g minus the (— logL) value calculated at the Arris where 
the (— logL) is at its minimum. It uses the values of A^ and Am^; determined above 
as inputs to the fit, and assumes a B^ fraction of /^s = 12%. The data prefer high 
values of Amg, with a favored value of 8 ps^^. The difference in likelihood for higher 
values of Ams is insufficient to exclude them, therefore a lower limit is set on Ams. 
Superimposed on the data is a 95% confidence level lower limit curve calculated using 
a 'fast' Monte Carlo. 

In constructing the limit curve, the likelihood differences, AlogL, for the fast 
Monte Carlo are calculated for 300 samples at various input values of Amg (2.0, 
4.0, 5.0, 6.0 and 7.0 ps~^), each with sample size equal to that of the data. If the 
Ams value is close to the point where the limit is set, 600 samples are used. The 
95% confidence limit is determined by locating the point below which lie 95% of the 
A logL values, calculated at the input value of Amg. The 95% confidence limit curve 
is then drawn through the points at different input Amg, as shown in Fig. |^(a). The 
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Fig. 12. Limit on Arris from (a) DELPHI and (b) OPAL Lepton-Jet charge method. 

data A log!/ curve intersects the limit curve at Arris = 6.1 ps~^. This point is taken 
as the 95% confidence level lower limit. The lower plot of Fig. ^(b) shows the result 
of performing this complete analysis with several different values of Jbs as discussed 
later in this section. 

It is important to check that there is indeed sensitivity at Arris = 6 ps~^. For 
this purpose, 800 Monte Carlo samples were generated at this value of Arris with the 
statistics of each sample again matching those of the ALEPH data. For each of these 
800 samples, the value of Arris at the minimum of the (— logL) curve is determined. 
The distribution of these minima are shown in Fig. |I0[ The figure clearly shows 
that in the majority of cases, the method does find the correct minimum, which 
demonstrates that there is indeed sensitivity at Arris = 6 ps^^. 

The data curve of Fig. ||(a) corresponds to an input of 12%. Fig. ^ shows 
the corresponding data curves obtained from the ALEPH data for various assumed 
Bs fractions. The figure demonstrates that the sensitivity to Arris increases as fss 
increases. Fig. P(b) shows the results of 95% confidence level lower limit in Arris as 
a function of fss- The limit varies from Anig > 5.2 ps"'^ at fss = 8% to Arris > 6.5 
ps~^ at /bs = 16%. 

The corresponding preliminary results using the Lepton-Jet charge method from 
Delphi! and OPAlI are shown in Fig |l|. The result from DELPHI is Ams > 
4.2 ps-i at 95% Confidence Level for Jbs = (10 ± 3)%, while from OPAL, the result 
is Arris > 3.3 ps~^ at 95% Confidence Level for fss = (12 ± 3.6)%. Taking the 
data curve in Fig. |l2|(b) literally, OPAL also excludes ranges of Arris between 6.3 and 
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Fig. 13. Limit on Anig from (a) ALEPH and (b) OPAL Lepton-Lepton method. 



7.9 ps ^ and above 19.6 ps ^ at 95% Confidence Level. At 97% Confidence Level, 
however, these exclusions disappear, so their significance is marginal. 



4-2. Lepton-Lepton Method 

This method has been used by the ALEPH0 and OPALi collaborations at LEP, 
and its schematic is that of Fig. |^(a) with the and replaced by Bs and Bg. As 
with the Bd analysis, the B decay vertex is formed by the secondary vertex including a 
high pt lepton, and the decay and production flavors are tagged by the signs of leptons 
on the flight distance and opposite sides. The results from the ALEPH and OPAL 

n 

collaborations are shown in Fig. |TB|. The preliminary result from ALEPHQ is Amg > 
5.6 ps^^ at 95% Confidence Level for = (12.2 ± 3.2)% while the published result 
from OPAlJ is Am, > 2.2 ps-^ at 95% Confidence Level for fe, = (12.0 ± 3.6)%. 



4-3. Lepton-Kaon Correlations Method 

This method has been used by the ALEPhEI collaboration, and its schematic 
is given in Fig. |T^(a). In order to enrich the sample with B, events, the analysis 
requires that a charged kaon from the primary vertex be identified. This kaon must 
have the opposite sign of the lepton or jet charge in the opposite hemisphere, in order 
to improve the tag of the production flavor. To enrich the decay vertex with Dg, the 
"charm vertex" is required to contain either zero or two kaons, or one kaon with a 
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Fig. 14. Limit on Amg from ALEPH with Lepton-Kaon Correlations method. 



charge opposite to the lepton. The decay flavor is tagged by the sign of the lepton 
from the decaying B meson, as in other methods. 

This selection yields 4436 lepton-kaon correlations, and enriches the Bs sample by 
a factor of 1.35. The method has a high tag rate of about 80%. The preliminary result 
for this measurement is shown in Fig. [l5(b) giving Arris > 4.0 ps^^ at 95% Confidence 
Level for fe^ = (12 ± 3)%. 



4.4- Estimation of B^ fraction in b events. 

There are currently two methods for determining the fraction of Bs in an inclusive 
lepton sample. The first method is from Dsi correlations. ALEPH has measured the 
product branching ratioEl: 

fBs ■ Br(5° ^ D-tvX) = 0.82 ± Om{stat)toil{syst)% (12) 



and from this, derivediaEl /^^ = 11.0 ± 2.8%. 

The second method uses the average time-integrated mixing parameter x = 
fssXs + fs^Xd- Assuming /^^ + fg^ + fss + fK, = 1, /b„ = fs^, and Xs = 0.5, 
the -Bs fraction is given by: 

f 2X - (1 - fAt,)Xd 

IBs : ■ l-'-'^j 

1 - 

UsingEI tb^ = 1.570 ± 0.049 ps, Am^ = 0.458 ± 0.020 from the LEP average as given 
in Fig. ^ and the T(4s) averaged Xd(^(4s)) = 0.167 ± 0.025, the world average of 
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Table 1. Summary of measurements of x- 



Measurement 


X 


LEP+SLDE 
(dileptons] 




0.1145 ± 0.0061 


CDF (e^)E 


r 


0.118 ± 0.008 ± 0.020 


CDF (/i^J 


i- 


0.118 ± 0.021 ± 0.026 


DO {fifiP 




0.09 ± 0.04 ± 0.03 


World Avera 


ge 


0.115 ± 0.006 



the time-integrated Bd oscillation parameter, Xd is calculated to be 0.170 ± 0.011. 
Table m gives the average mixing parameter, x = 0.115 ± 0.006. The baryon fraction, 
/aj, is derived from Ac-lepton and A-lepton correlations, in analogy with Eq. ([T^) 
above. Using! Br(Ac AX) = 35 ± 11%, the LEP measurementJH3 can be averaged 
to yield 

fA, ■ Br{Ab AcXiu) = 1.67 ± 0.30%, (14) 

with common systematic effects taken into account. Following the method given in 
Ref. 20, the baryon fraction is then calculated to be = 12.8 ± 3.9%. 

With these inputs, Jbs from the second method is then 9.9 ± 1.9%. An average 
of the two methods then yields a final estimate of the fraction of Bs mesons produced 
in Z ^bb decay, Jbs = 10.2 ± 1.6%. 

4-5. Summary of lower limit for ArUg 

Table ^ summarizes the lower limits at 95% Confidence Level placed on Arrig from 
the LEP experiments. These limits on Arris are computed using different techniques, 
and there is currently no combined result which takes correlated statistical and sys- 
tematic errors into account. Thus, the best limits on Bg oscillation. Arris > 6.1 ps~^ 
for fss = 12% and ArUg > 5.6 ps~^ for /^g = 10% using the Lepton-Jet charge 
method by ALEPH are taken as the current limits on Arris. Defining Xg = ArrisTBs 
whereill tbs = 1-58 ± 0.10 ps~^, the values of Xs are shown in Table 3 by shifting 
the central value of down by la. Using the world average central valuesEl of the 
quantities in Eq. (|^) and including their uncertainties by shifting the values by la to 
the conservative side, yields the ratios Arris/ Am d and |^s/Vtd| as shown in Table 3. 
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Table 2. Summary of Limits on Atos at 95% C.L. 





s V 1-^^ y 


fn 

J i>s 


ALEPH (91-94) 


> 6.1 


12% 


flent/O i) 


> 5.6 


10% 


ALEPH (91-94) 


> 5.6 


12±3% 


(lept/lept) 






ALEPH (91-94) 


> 4.0 


12±3% 


(lept/K+O t) 






DELPHI (91-94) 


> 4.2 


10±3% 


(lept/Qj) 






DELPHI 


> 1.5 




(DJ/Qj) 






OPAL (91-94) 


> 3.3 


12.0±3.6% 


(lept/Qj) 






OPAL (91-93) 


> 2.2 


12.0±3.6% 


(lept/lept) 







Table 3. Constraints on physical quantities resulting from measurements of Am^ and Arris- 





fBs = 12% 


fBs = 10% 


Arris 


> 6.1 ps~i 


> 5.6 ps^^ 


Xs 


> 9.0 


> 8.3 




> 12.8 


> 11.8 


\Vts/Vtd\ 


> 2.8 


> 2.7 



5. Conclusion 

In summary, by studying the time-dependence of B^-B^ oscillations, recent ex- 
periments have given 

(i) an accurate value for Am^, the mass difference between {B(i)l and (-8^)5; and 

(ii) a lower bound for Arris, the mass difference between (-Bs)l and {Bs)s- 

The results from the ALEPH, DELPHI and OPAL Collaborations at LEP of 
CERN and the CDF Collaboration at the Tevatron Collider of Fermilab are summa- 
rized in Table ^. In particular, 

^^ = 85.8 ±3.6. (15) 
ArriK 

The impact of the results presented here is shown in Fig. 15 in the (Am^, Arris) 
plane together with the region allowed by the Standard Model0. 
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Table 4. Summary of Am^ and Arris results. 



Mass Differences for the Long and Short Eigenstates 




Am (ps ^) 


Am (eV) 


ArriK 


(5.33 ±0.03) X 10-=^ 


(3.51 ±0.02) X 10-6 


Arrid 


0.457 ±0.019 


(3.01 ±0.13) X 10-^ 


Am, (95% C.L.) 






fBs = 12% 


> 6.1 


> 4.0 X 10-3 


/Bs = 10% 


> 5.6 


> 3.7 X 10-3 




Fig. 15. Constraints on the (Ato^, Arrrg) plane. 
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